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Abstract 

We present an analytical fit to the growth function of the dark matter perturbations 
when dark energy perturbations are present. The growth index 7 depends upon the 
dark energy equation of state w, the speed of sound of the dark energy fluctuations, 
the dark matter abundance and the observed comoving scale. The growth index 
changes by (9(5%) for small speed of sound and large deviations of w from —1 with 
respect to its value in the limit of no dark energy perturbations. 
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1 Introduction 

Current observations of Type la supernovae luminosity distances indicate that our 
Universe is in a phase of accelerated expansion [1]. Various proposals have been 
put forward to explain the present acceleration of the Universe. One can roughly 
distinguish two classes. On the one hand, the acceleration might be caused by the 
presence of dark energy, a fluid with negative equation of state w. This may be pro- 
vided by a tiny cosmo logical constant which is characterized by w = —1 or by some 
ultralight scalar field whose potential is presently dominating the energy density 
of the Universe. This is usually dubbed quintessence [2] (see [3] for a comprehen- 
sive review). On the other hand, the acceleration might be due to a modification 
of standard gravity at large distances. This happens in f(R) theories [4] and in 
extra- dimension inspired models, like DGP [5]. Understanding which class of mod- 
els Nature has chosen will represent not only a breakthrough in cosmology, but 
also in the field of high energy physics. 

Mapping the expansion of cosmic scales and the growth of large scale structure 
in tandem can provide insights to distinguish between the two possible origins 
of the present acceleration. For such reason, there has been increasing interest 
in analysing the time evolution of the dark matter perturbation. Several recent 
works deal with characterizing the growth of dark matter perturbations in different 
frameworks [6-20]. 

The evolution of the growth function of dark matter perturbations g = 5 c /a , 
which is the ratio between the perturbation 5 C and the scale factor of the Universe 
a, can be parameterized in a useful way using the growth index 7 [21], defined in 
Eq. [19j In a pure matter-dominated Universe, g does not evolve in time (remains 
equal to one) and 7 is zero. However, in the presence of a dark energy background, 
g changes in time, 7 is different from zero and its value can be approximated by 

7 = 0.55 + 0.05 [1 + w(z = 1)], (1) 

which provides a fit to the evolution of g to better than 0.2% for — 1 < w and a 
broad range of initial conditions for the dark matter abundance [21]. Typically, 
the growth index in modified gravity models turns out to be significantly different 
(for instance 7 ~ 0.68 for DGP [21]) and therefore it is in principle distinguishable 
from the one predicted for dark energy modelgj. The available data on the growth 
of structures are still poor and there is a long way to go before we can talk about 
precision cosmology in this respect. The methods developed to study the growth of 
structure involve baryon acoustic oscillations, weak lensing, observations of X-ray 
luminous clusters, large scale galaxy surveys, Lyman-a power spectra and the 
integrated Sachs- Wolfe effect on the Cosmic Microwave Background. There are 
however various works that use these kind of techniques to place constraints on the 
growth index (and some also on the equation of state of dark energy) as well as 



^However, see [9] and [20]. 



forecasts for its determination based on future observations [22-31]. In particular, 
it is found in [30] using Bayesian methods that a next generation weak lensing 
survey like DUNE [32] can strongly distinguish between two values of 7 that differ 
by approximately 0.05 . The authors of [23] made a forecast for the same kind 
of satellite proposal and concluded that it will be possible to measure the growth 
index with an absolute error of about 0.04 at 68% confidence level. In [24] a slightly 
bigger error of 0.06 at the same confidence level is given for a forecast based on 
baryon acoustic oscillations. Finally, for a combination of weak lensing, supernovae 
and Cosmic Microwave Background data an error of about 0.04 is estimated in [22] 
after marginalizing over the other cosmological parameters. Since the growth index 
is approximately equal to 0.55, the nearest future observations should be able to 
determine it with a relative error of around 8%. 

While much effort has been put into determining the value of the growth index 
in dark energy and in modified gravity models, less attention has been devoted 
to the possible effect on 7 of non-vanishing dark energy perturbations. The lat- 
ter do not affect the background evolution, but are fundamental in determining 
the dark energy clustering properties. They will have an effect on the evolution 
of fluctuations in the matter distribution and, consequently, on 7. While mini- 
mally coupled scalar field (quintessence) models commonly have a non-adiabatic 
speed of sound close or equal to unity, and therefore dark energy perturbations 
can be neglected for them; other non-minimal models, for instance the adiabatic 
Chaplygin gas model, motivated by a rolling tachyon [33], have a speed of sound 
which is approximately zero. Observational implications of dark energy perturba- 
tions with a small speed of sound in a variety of dark energy models have been 
recently discussed in k-essence [34,35], condensation of dark matter [36] and the 
Chaplygin gas, in terms of the matter power spectrum [37, 38] and combined full 
CMB and large scale structure measurements [39,40]. Let us also emphasize that 
dark energy perturbations may not be consistently set to zero in perturbation the- 
ory [41] even if w = — 1. Indeed, it is unavoidable that dark energy perturbations 
are generated, even if set to zero on some initial hypersurface, due to the presence 
of a non-vanishing gravitational potential. Therefore, the expression (TjQ) rigorously 
holds only in the physical limit in which the speed of sound is very close to unity 
(if w 7^ —1) so that dark energy perturbations are sufficiently suppressed. 

In this Letter we study the effect of dark energy perturbations on the growth 
index 7. Our main motivation is to understand if the introduction of the new 
degrees of freedom introduced by dark energy perturbations imply changes in 7 
large compared to the forecasted errors A7 ~ O(0.04) (at 68% confidence level). 
Following the common lore, see for instance [39], and to simplify the analysis, we 
will assume that the speed of sound associated with the dark energy perturbations 
and the equation of state do not change appreciably in the proper time range and 
that the dark energy perturbations have no shear. This is a good approximation in 
linear perturbation theory for dark energy models with a scalar field. Under these 



assumptions, we provide an analytical formula for the growth index 7 as a function 
of the speed of sound, the equation of state w, the dark matter abundance and the 
comoving scale. As we will see, in the presence of dark energy perturbations, the 
growth index differs from the corresponding value without dark energy perturba- 
tions by an amount which is comparable to the realistic forecasted errors, especially 
for small speed of sound and w significantly different from — 1. This opens up the 
possibility that the presence of dark energy perturbations may leave a significant 
imprint on the growth function of dark matter perturbations. 

The Letter is organized as follows. In Section 2 we summarize our framework 
and provide the necessary equations for the perturbations at the linear level. In 
Section 3 we discuss the growth index and in Section 4 we give our results and 
summarize of our work. 



2 The basic equations 

In this section we shortly describe how to obtain the second order differential equa- 
tions describing the evolution of the coupled linear perturbations of dark matter 
and dark energy in a spatially flat Friedmann-Lemaitre- Robertson- Walker (FLRW) 
background. We will closely follow [42] and [39] and work in the synchronous gauge 
for convenience. With this choice the perturbed metric in comoving coordinates 
reads 

ds 2 = a 2 (r) [-dr 2 + (<% + h tj ) dx l dx j ] , (2) 

where hij encodes the perturbation and can be decomposed into a trace part h = h\ 
and a traceless one. The background equations are simply 

3H 2 = 8nGa 2 p, (3) 

2H' = -n 2 (1 + 3 W n x ) , (4) 

where G denotes Newton's constant, p = p c + p x is the total energy density, the 
comoving Hubble parameter is 7i = a' /a, primes denote derivatives with respect to 
the comoving time r and we define the time varying relative dark energy density 
as Q x = Px/p- The bars indicate homogeneous background quantities and the 
subindexes ' c ' and V refer to dark matter and dark energy respectively. We assume 
that the equation of state of dark energy, w, is a constant and that the dark energy 
and the dark matter do not interact. The divergence of the dark matter velocity 
in its own rest frame is zero by definition and therefore in Fourier space we have 

S' c + lh' = 0, (5) 

where 

Sp c = p c S c , (6) 



is the energy density perturbation of dark matter. The speed of sound of a fluid 
can be defined as the ratio [39] 

c* = - (7) 

where we have introduced 5P, the pressure perturbation of the fluid. It is important 
to recall that the speed of sound defined in this way is a gauge dependent quantity. 
However, the speed of sound is gauge invariant when measured in the rest frame 
of the fluid. The pressure perturbation of a dark energy component with constant 
equation of state can be written in any reference frame in terms of its rest frame 
speed of sound c s as follows 

5P X = c 2 s 6p x + 3H(1 + w) (c 2 s - w) pj^ , (8) 

where 9 X is the dark energy velocity perturbation and k the inverse distance scale 
coming from the Fourier transformation. Then, taking into account this expression 
and the relation 

h" + Hti = 8irGa 2 (5T ° - 5T[) , (9) 

where T^ is the energy-momentum tensor, one can differentiate ([5]) with respect to 
r and make use of the background evolution ([3]), (@J to find the equation for the 
dark matter energy density perturbation [39] 



^ + m' c -ln 2 n c 5 c = ln 2 n x 

2 2 



(1 + 3c 2 ) 5 x + 9(l + w)H (c 2 - w) || 



(10) 



The time derivative of the dark energy density perturbation in the dark matter 
rest frame is [39] 

6* = -(1 + «0 { [k 2 + 9 (c 2 s - w) H 2 ] |f - A - 3H(c 2 - w)5 x (11) 

and the time derivative of the divergence of the dark energy velocity perturbation 
in the case of no anisotropic stress perturbation is 

§ = -(i-^)4 + T iU. (12) 

Differentiating ( TTTjl with respect to the comoving time and combining ( fill) and ( 1121) 
with the background equations into the resulting expression one gets 

b" £ + [3(c 2 s -w)H-F]6' x 

+ U 2 s k 2 -^(c 2 s -w)H [(1 + 3wQ x - 6c') U + 2T] 1 S x 

= {l + w)%-(l + w)F6' e , (13) 



where 



c 2 — w 



y = -9 (1 + 3*>.) k2 + 9 '., _ w) H2 H* - (1 - UQH . (14) 
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Equations ffTUl) . ffTTl) . f|T3|) and (fT4"l) allow us to describe the evolution of linear 
perturbations of dark matter and dark energy as functions of time in a FLRW 
background. Initial conditions are given at the redshift z mr = 3200, which approx- 
imately corresponds to the time of matter-radiation equality. Since we consider 
non-interacting fluids to describe the dark matter and dark energy, their energy 
densities satisfy: 

p c ' + 3Hp c = 0, (15) 

pj + 3(1 + w)Hp x = 0. (16) 

We choose adiabatic initial conditions 

$x(mr) = (1 + w)5 c{mr ) . (17) 

Furthermore, we assume zero initial time derivatives of the matter and dark energy 
perturbations. This is consistent with the fact that at early times (both in the 
radiation and matter dominated periods) the equations of the perturbations admit 
the solution S x oc (1 + w)5 c oc r 2 [39] as can be checked with ( ITU]) . ( TTTT) . ( 1T31) and 
( TT4")) and and we have set the initial conformal time to zero. In fact we can even use 
non-zero initial velocities and consider non-adiabatic initial conditions; our results 
on the growth index are robust under these modifications. For the background we 
consider the present (i.e. at z = 0, a = 1) value of the relative energy density of 
dark matter in the range (0.25,0.30) and f2° = 1 — f2°. In our computations we 
do not include a specific baryon component. We have checked that the effect of 
adding baryons on the growth index can be at most as big as 0.2%, which is much 
smaller than the 8% accuracy forecasted for the near future experiments. 



3 The growth index 

The growth of matter perturbations has been studied neglecting the effect of dark 
energy perturbations through the behaviour of the growth function [43] 

g = - (18) 

a 

as a function of the natural logarithm of the scale factor. It is possible to fit g 
using a simple parameterization that defines the growth index 7 and depends on 
the relative energy density of dark matter 

g(a) = gfa) exp / (fi c (a) 7 - 1) 4? • (19) 

The growth function depends on the scale k, the sound speed c 2 s and the equation of 
state w. This dependence is embedded in the growth index 7 which therefore from 
now on has to be understood as a function of these parameters. The growth factor 



g can be normalized to unity at some a% > ai mr \ deep in the matter dominated 
epoch where S c ~ a. The growth index 7 is normally taken to be a (model- 
dependent) number whose best fitting value for standard gravity and no dark energy 
perturbations is around 0.55, see equation flTJ). This result is obtained from the 
equation 

% + H8' e -~'H?n e 6 e = 0, (20) 

with no dark energy perturbations, instead of the system of second order differential 
equations that includes S x . 



It is important to remark that it is not possible to reduce the system f fTOj) . (fTTj) . 
( ITS]) and (fll|) to (f2"U|) by setting S x = or with any particular choice of the parame- 
ters. Those equations show that even if the dark energy perturbation is set to zero 
initially it will be generated at later times. The effect of dark energy perturbations 
should be included in the analysis of the growth history for consistency. The growth 
of dark matter perturbations depends not only on w (which already enters in (1201 
through Q c and 7i) but also on the other two parameters appearing explicitly in 
the differential equations that control the evolution of the perturbations, i.e. k 
and c 2 s . The reason for the dependence of the dark matter perturbations on the 
sound speed of dark energy is clear from the previous discussion and the definition 
(j7p. In contrast to equation (1201 . the dependence on the comoving momentum now 
appears explicitly as an effect of a non-vanishing speed of sound. 

Given the numerical solution for the dark matter perturbation evolution, the 
definition ( fl9l of the growth index can be used to compute 7 exactly: 

T^Mit)- (21) 



In the next section we will use this equation together with ( fTOl . ( TTTT) . ( fT3l) and ( flip 
for obtaining our results. Obviously 7 will be a function of a and it will depend on 
k, ct, w and Q®. as well. 



s' 



In our analysis we consider w in the reasonably broad range (—1, —0.7). We 
choose not to allow the possibility that the equation of state of dark energy can be 
smaller than —1. As for k, the values of interest are the ones for which there is large 
scale structure data on the matter power spectrum [44]. This goes approximately 
from O.Olh Mpc -1 to 0.2h Mpc -1 , including the nonlinear part of the spectrum 
which becomes so at roughly 0.09h Mpc -1 . The scale that corresponds to the 
Hubble size today is 2.4 10 -4 Mpc -1 and if we normalize it to 7io = 1 , the range of 
k we will focus on (discarding the nonlinear part of the spectrum) is approximately 
(30, 270) in units of 7i • Notice that the lower k value roughly gives the position of 
the baryon acoustic oscillation peak that can be used for constraining the growth 
index [24]. Finally, regarding the sound speed of dark energy, we restrict c 2 s to 
be positive and smaller or equal than unity as currently the bound is very weak 
[39,45-49]. 
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Figure 1: j(z = 1) as a function of w is shown for four values of c s . Red curves 
correspond to k = 0.050 /iMpc -1 and blue dashed ones to k = 0.078 /iMpc -1 . 



4 Results and discussion 



In this section we present a combination of numerical results and an analytical for- 
mula for the growth index 7 as a function of the relevant cosmological parameters. 

In Figure [I] we plot the growth index at z = 1 versus w for several values of 
the speed of sound of dark energy and two different scales. Notice that the curves 
for the two different values of the comoving momenta coincide for c s — 1 and in 
the limit of very small speed of sound. The figure indicates that the dark energy 
speed of sound and the scale determine whether 7 grows or decreases as a function 
w at a given redshift. This is one of the reasons why having a more complete 
parameterization than (CQ) is important. Choosing another redshift would have the 
effect of an overall shift of the merging point at w = — 1 together with modifications 
in the curvatures of the lines. 



To gain some insight on the change of the value of 7 from c 2 = 1 to c 2 <^C 1 , we 
observe that, in the limit (? s ~ and from Equation (!T2l) . the dark energy velocity 
perturbation promptly decays in time. One is left with the following solution for 

) +(H-1«)0 

where the dot stands for differentiation with respect to In a. As a first approx- 
imation, we can solve Equation fl22l) plugging in the dark matter perturbation 
S c (Cg = 1) obtained taking c 2 s — 1, which for this purpose corresponds to the case 



3-u; 



~- 3 „-l 



S c (c 2 s = l) da, 
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Figure 2: Relative error as a function of w between the exact numerical result for 
7(2 = 1) with very small dark energy speed of sound and the approximation 7 ap at the 
same redshift based on Equations (|22p and (|10p with zero 8 X . The figure has been done 
10- 6 , n° c = 0.30 and k = 0.050 /iMpc" 1 . 



for c 2 s 



in which no dark energy perturbations are present. From Equation fllOl) . it is 
clear that the dark energy perturbations provide an extra source for the dark mat- 
ter pertrubation growth. We then solve numerically Equation (fTUj) with this new 
known source and 9 X = . The difference between the true value of 7 and the one 
obtained with such an approximation is plotted in Figure [2j 

In Figure [3] we show the growth index at z — 1 versus log 10 c s for different 
values of the equation of state of dark energy and two scales k. From this plot it is 
clear that the effect of changing the scale is an overall shift along the log 10 c s axis. 
Notice that the intersecting points for the two sets of lines have the same value of 
the growth index, 7 ~ 0.547, which corresponds to the merging point in Figure [H 

The redshift dependence of the growth index has already been studied without 
taking into account dark energy perturbations [50] concluding that d'y/dz ~ —0.02 
at z = ; being this value nearly independent of z for a given f2° . However, 
including dark energy perturbations, we find that it is actually the derivative of 
7 with respect to the scale factor a which is constant. Therefore the redshift 
dependence of the growth index can be better modeled with a \jz term plus a 
constant term. We will later see that the growth index actually has an almost 
constant slope as a function of the scale factor when dark energy perturbations are 
taken into account. 
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Figure 3: 7(2 = 1) as a function of log 10 c s is shown for four values of w. Red curves 
correspond to k = 0.03 /iMpc~ and blue dashed ones to k = 0.08/iMpc~ . 



Our next step is to obtain an analytical parameterization of the growth index 
as a function of the cosmological parameters. We start with the following generic 
ansatz: 

-f(fl° c ,c s ,k,w,a) = 7^ (Q° c ,c s ,k,w) +( (Q° c ,c s ,k,w) [a - a eq (0°,u>)] , (23) 



where a eq is the value of the scale factor at which "dark equality" {Vt c 
takes place: 

' 1 



"eq 



n°. 



-1 



n x = 1/2) 



(24) 



We want to fit the growth index for a in the interval [a eq , 1] which approximately 
corresponds to a redshift z G [0, 0.55] for the ranges of the equation of state of dark 
energy and its relative energy density that we consider. Ideally one would wish 
to be able to use (|2T|) and the equations for the perturbations to infer completely 
the analytical dependence of •y eq (fi°, c s , k, w) and ( (Q°, c s , k, w) in their variables. 
This turns out to be difficult and we find it efficient to make a numerical fit directly. 
The generic form (|23|) . which can be viewed as a first order Taylor expansion in 
the scale factor, is motivated by the nearly zero variation of d7/da. The choice of 
a eq as the point around which we make the expansion is a convenient one, but the 
fit could in principle be done taking a model independent value of a as the fiducial 
point. We use the same ansatz to fit 7^ and 70 , which is the growth index at 
a = 1, and doing so we directly obtain the slope £ from ( |23l : 



C(n° c ,c s ,k : w) 



70 - leq 



(25) 
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Figure 4: ^ eq (c s , k, w) versus log 10 c s for different combinations of the pair {k, w} : 
A = {0.08 h Mpc -1 , -0.95} , B = {0.02 h Mpc" 1 , -0.7} , C = {0.04 h Mpc" 1 , -0.87} and 
D = {0.06 h Mpc -1 , —0.75}. Red lines are the exact numerical result and blue dashed 
ones the corresponding fits. 



In particular, we assume the following parameterization for 7 eg and 70 : 



7j (c s , k, w) = hj(w) tanh 



(to* 



10 c s 



<?#)) 



Tj(w) 

hj(w) 



+ fj(w) , j = {eq,0} . (26) 



Notice that we have taken ~/ eq and 70 to be independent of Q® and we incorporate 
this assumption in our notation, so we will refer to jj (c s , k,w) from now on. The 
functions fj(w), gj{k) , hj(w) and rj(w) are polynomials in their variables. It 
turns out that the fit obtained with this procedure can be importantly improved 
with the addition of a polynomial correction to ( that depends on fi°, so finally: 

j(fl° c ,c s ,k,w,a) = 'j eq (c s ,k,w) 

+ [C (n° c , c s , k,w)+ V (0°)] [a - a eq (fi°, u;)] . (27) 

The set of equations (1251) . (1261) and (1271) constitute the full fitting formula for 
the growth index. The resulting nine polynomials through which the fit can be 
expressed are the following: 

f eq { w ) = 4.498 -10 -1 -2.176 • 10 -1 w - 1.041 • 10 -1 w 2 + 5.287 • 10~ 2 w 3 

+ 4.030 • 10~ 2 w A , (28) 

f (w) = 4.264 -10 -1 -3.217 ■ W' 1 w - 2.581 • 10 -1 w 2 - 5.512 • 10" 2 w 3 

+ 1.054 • 10~ 2 w A , (29) 
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9eq (k) = -5.879- KT 1 -2.296 • 1(T 2 k + 2.125 • 10~ 4 k 2 -1.177 • 1(T 6 A; 3 

+ 3.357 • 1(T 10 k A - 3.801 • l(T 12 /c 5 , (30) 

go (k) = -6.401 • 10 _1 -2.291 -10~ 2 k + 2.119 • 10 -4 k 2 -1.173 • 10~ 6 k 3 

+ 3.344 • 10~ 10 k A - 3.787 • 10~ 12 £; 5 , (31) 

h eq (w) = 1.759- 10 _1 + 4.066- 10 _1 iu + 3.254 • 10" 1 w 2 + 9.470 • 10~ 2 u> 3 , (32) 
h eq (w) = 2.008 • 10" 1 + 4.644 • 10" 1 w + 3.713 • 10" 1 w 2 + 1.076 • 10" 1 w 3 , (33) 

r eq (w) = 5.158 • KT 1 + 1.203u> + 9.697 • 10" 1 w 2 + 2.827 • 10" 1 w 3 , (34) 

r ( w ) = 6.093 • lO -1 + 1.435u> + 1.1668w 2 + 3.412 • 10" 1 w 3 , (35) 

ri(n°) = 8.037 • 10~ 3 + 4.676 • 10~ 2 0° - 2.829 • 10" 1 (Q° c ) 2 . (36) 

The truncation of the coefficients above has been done in such a way that the figures 
in the Letter can be reproduced and that the maximum relative error between 
the numerical value of 7 and the fitting formula does not exceed 0.2% for any 
combination of the parameters. In fact, this error turns out to be much smaller for 
generic choices of the parameters. 

In Figured] we show 7 eg (c s , k, w) versus the decimal logarithm of c s for several 
combinations of k and w. The red curves represent the exact numerical growth 
index and the blue dashed lines are the corresponding fits. In Figures [51 [6] and [7] 
we show 7 (fi°, c s , k, w, a) versus the scale factor for several values of w, f2° and 
k respectively, as explained in the captions. The other parameters are kept fixed. 
The colour code, as in Figure HI is that the red curves represent the exact numerical 
growth index and the blue dashed lines are the corresponding fits. These figures 
are meant to illustrate the goodness of fit for several choices of the parameters. 

Equations (25)- (36) offer an analytical expression for the growth index in terms 
of the relevant cosmological parameters in the case in which dark energy pertur- 
bations are present. The case without dark energy perturbations is reproduced by 
assuming c 2 = 1 . The analytical parameterization fits the numerical results in the 
assumed range of parameters to a precision of 0.2% (in the worst cases) or better 
for the growth index. Our findings show that 7 can vary from 0.55 by an amount 
A7 as large as ~ 0.03 . We have checked that this result holds for any redshift 
between z eq (at the time of dark equality) and z — 1. This difference is of the same 
order of magnitude of the 68% c.l. forecasted error band. The predicted value of 7 
may differ by this amount from the value without dark energy perturbations if the 
speed of sound is tiny and if the equation of state substantially deviates from — 1. 
This opens up the possibility that a detailed future measurement of the growth 
factor might help in revealing the presence of dark energy perturbations. Finally, 
let us reiterate that our results have been obtained under the assumption that c 2 
and w do not evolve in time, at least for mild values of redshift. Furthermore, we 
have assumed that the dark energy perturbations have no shear. 
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Figure 5: 7 (0°,c s , k,w,a) versus a for k = 0.033 /iMpc" 1 , 0° = 0.27 and c 2 s = 0.01 . 
Different values of w are chosen as shown in the figure. The red lines are the numerical 
results from the differential equations and the blue dashed ones are the fits to them. 
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Figure 6: 7 (fi°, c s , k, w, a) versus a for k = 0.03 hMpc' 1 , w = -0.92 and c 2 = 0.0036 . 
The value of Q® runs between 0.25 and 0.30 in steps of 0.01 from top to bottom of the 
figure. The red lines are the numerical results from the differential equations and the 
blue dashed ones are the fits to them. 
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Figure 7: 7 (fi°, c s , k, w, a) versus a for w = —0.80 , c 2 s = 0.01 and 0° = 0.27 . The scale 
k in units of h Mpc -1 takes the values {0.023, 0.027, 0.037, 0.067} from bottom to top of 
the figure. The red lines are the numerical results from the differential equations and the 
blue dashed ones are the fits to them. 
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